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Abstract. We consider a fully-quantized model of spontaneous emission, 
scattering, and absorption, and study propagation of a single photon from 
an emitting atom to a detector atom both with and without an intervening 
scatterer. We find an exact quantum analog to the classical complex analytic 
signal describing an electromagnetic wave scattered by a medium of charged 
oscillators. This quantum signal exhibits classical phase delays. We define a time 
of detection which, in the appropriate limits, exactly matches the predictions 
of a classically defined delay for light propagating through a medium of charged 
oscillators. The fully quantized model provides a simple, unambiguous, and causal 
interpretation of delays that seemingly imply speeds greater than c in the region 
of anomalous dispersion. 



PACS numbers: 32.80.-t, 42.25.Bs, 42.50.Ct 

1. Introduction 

The subtleties of appropriately defining the speed of light pulses in dispersive media 
have been investigated for over a century. Although a careful classical analysis was 
performed early in the twentieth century (see, for example, the papers collected in 
reference P), recent observations of anomalously slow and anomalously fast 
speeds of light have returned studies of light propagation to the pages of contemporary 
physics journals, and a recent review so-called superluminal propagation is contained 
in 1^1 . In this paper we study the scattering of single spontaneously emitted photons 
and develop a simple fully-quantized microscopic model of light propagation in a 
dielectric medium. We also compare this quantized model to a simple microscopic 
classical model of scattering from charged oscillators. Although our models are limited 
to the simplest case of weak scattering in linear dielectrics, they serve as a basis for 
understanding propagation in more complex media at the quantum level. 

Classical pulses can be understood as the superposition of the incident field with 
fields re-radiated by the atoms of the medium. This point of view has been articulated 
clearly by Feynman 0, with recent elaborations by Sherwood |S] and Milonni S'. 
The re-radiated fields are phase-shifted and attenuated relative to the incident field, 
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and the leading edges of all fields propagate at the vacuum speed of light c. We 
identify quantum mechanical quantities that have a striking quantitative parallel in 
the complex analytic signal describing fields in the classical model. Although the state 
vector describing the system has no absolute overall phase, the scattering-induced 
phase shift of the classical field is manifested in the description of the quantum field 
as well as in the probability amplitude for atoms subsequently excited by the field. 
The pulse re-shaping of classical fields that leads to group delays is also evident in 
our quantum model, and results in delays that are equivalent in the classical and 
quantum models. All delays are clearly due to the superposition of quantum effects 
that propagate at the vacuum speed of light c, in spite of the appearance of delays 
that seemingly imply so-called superluminal velocities. The parallels between classical 
and quantum scattered fields facilitate the understanding of delays of single photons 
in terms of classical concepts. 

The classical medium we consider consists of charged point particles attached 
to their equilibrium positions by linear restoring forces. This model is considered in 
undergraduate textbooks [3 ^j. The classical medium has an index of refraction 
which for a dilute collection of scatterers depends linearly on the density of the 
scatterers, and the index of refraction determines the phase and group velocities of 
radiation in the medium. 

Our quantum system consists of a set of two-level atoms at fixed positions 
interacting with the quantized modes of a one-dimensional multimode optical cavity. 
An initially excited atom spontaneously emits radiation into the modes of the cavity 
(a photon); the radiation is scattered by a second atom; and a third atom serves as 
a detector. We find an analytical expression describing the excitation of the detector 
atom, and an expression for the expectation value of the space- and time-dependent 
electromagnetic field intensity. We then identify an exact quantum analog to the 
complex analytic signal of a classical field. The scattering-induced time delays in 
the quantum mechanical quantities thus have direct parallels in classical delays. We 
also find a simple, unambiguous, and causal interpretation at the quantum level of a 
superluminal velocity that occurs in the region of anomalous dispersion very near to 
an atomic resonance. 

We note that the results obtained analytically in this paper may also be observed 
using straightforward numerical techniques like those used in references jlllll2l[T^ll4j . 
Our work bears a particularly close relationship to the numerical work reported in I15j , 
in which quantum interference effects are observed in the scattering of single-photon 
wavepackets. 

We review the classical model in section|21 and in sectionOwe discuss our quantum 
mechanical model and present our analytical solution. (We leave the details of the 
derivation to the Appendix.) We highlight the similarity of the quantum and classical 
dynamics in section 0] and in section El we investigate a delay that seemingly implies 
superluminal propagation in both the classical and quantum models. 

2. Review of Scattering from Classical Oscillators 

A simple classical model of the interaction of radiation with matter consists of an 
electromagnetic plane wave driving electrons attached to molecules by linear springs 
[3 . This simple model is able to account for the attenuation and phase shift of 
waves transmitted through a dilute material comprised of such oscillators. 

For a plane wave normally incident on a dielectric slab with thickness Az and 
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real index of refraction n, the delay in the arrival of a point of constant phase on the 
far side of the medium (compared to a wave traveling in vacuum) is determined by 
the phase velocity = u/k = c/n. This delay is given by 

Az Az Az, , 
At^ = = (1) 

c c 

The delay in the arrival of the peak of a modulation envelope of a quasi-monochromatic 
pulse is determined by the group velocity Wg = duj/dk = c/{n + Lo^), and is given by 

Atfr^ = {n-l+ u— . (2) 

^ c c \ dco J 

For pulses that are not sufficiently monochromatic, the simple concepts of phase 
and group velocity are inadequate to characterize all of the effects of pulse-reshaping 
as the field propagates. Several other velocities and delays have been developed (see, 
for example, references I17| l and in this paper we focus on a delay determined by 
the "temporal center of gravity" of the field intensity of a pulse at a fixed position z 
"downstream" from the dielectric, i.e.. 

At,. ^ fjjn^i^) _ fjj£i^:^\ 

\ fS{z, t)2 dt J ^j^j^ dielectric V / '^(^' J no dielectric 
This is closely related to concepts used to define the centrovelocity in JS]. For quasi- 
monochromatic pulses far from resonance this delay is equivalent to the group delay, 
but in general it is necessary to calculate explicitly the field £ in order to determine 
Atg2. For the spontaneously emitted pulses with Lorentzian spectrums that are 
considered in this paper we will show that the "temporal-center-of-gravity" delay 
happens to be equal to twice the group delay. 

In order to find classical expressions for the index of refraction and absorption 
we follow the development outlined by Feynman We consider a monochromatic 
linearly-polarized plane wave of frequency u normally incident on an infinite thin slab 
of material with density of oscillators N and thickness Az. The plane of the slab is 
normal to the z direction, and located at z = 0. The oscillators have natural frequency 
(jJq, mass m, instantaneous speed and are assumed to experience a damping force 
proportional to the first derivative of acceleration, which for sinusoidal oscillations is 
equivalent to a damping proportional to —v. The field on the far side of the slab is the 
superposition of the incident field and the field scattered by the material. In the steady 
state the scattered field is proportional to the incident field, with an amplitude and 
phase shift given by simple resonance theory. The infinite slab geometry considered 
here creates an effectively one-dimensional model that matches the simple quantum 
mechanical system we introduce in the following section. 

If the incident field is Ei = Eq exp [— iw(t — z/c)], then the scattered field is [71111]) 
qNAz 
2eoc 

= i/ . 2 T ■ ^ -^i' (4) 

where in the last line we have introduced the dimensionless parameter / = 
N Az / {2meocj) which characterizes the magnitude of the scattering.. The total 
transmitted field Et on the far side of the material is 

Et = Ei + E, 

^LOq — LO^ — iju ^ 



= Ei 



1 + 



(5) 
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Figure 1. Phase and group delays of a plane wave propagating through a medium 
of classical oscillators. The group delay is negative in the region of anomalous 
dispersion, i.e., for —7/2 < (5 < 7/2. 



In comparing classical predictions with those from a quantum mechanical model 
we will be interested in the limit in which the detuning from resonance is small 
compared to the natural frequency of the oscillator, i.e., jcj — wqI ^ ^o- Taking 
advantage of this condition we rewrite the transmitted field of equation ([SJl in terms 
of the detuning from resonance 

J^w-wo, (6) 

and we approximate other occurrences of lo with wq, giving 



7 



^2(5 + i7 

For weak scattering, i.e., when / ^ 1, the transmitted field is approximately 



(7) 



1-/ 



7 



E\ exp 



(4^2 j_ ^,2^ 



■7 
2 



exp 



7(5 



(4(52+72) 

7(5 



(4(52 + 72) 

The phase-shift and attenuation given by equation 

7V(j2 7(5 

Tl ^ 1 — 

TOeoCt;o7 4^2 + 72 
and intensity absorption coefficient 



^i/ 



(4,52+ 72) _ • 
lead to the index of refraction 

(9) 



iVg2 



7 



meQCj (4:0'' + 7^] 

Using the expression for the index of refraction given in equation and 
assuming 5 ^ loq, the phase and group delays can be written 



2/ 



AtA ~ — - 



7(5 



uo \452 + Y 



(11) 
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Figure 2. Quantum mechanical model consisting of two-level atoms at fixed 
positions in a one-dimensional multimode optical cavity. 



(12) 

The phase and group delays predicted from the classical model are illustrated in 
figure^ The group delay is positive away from resonance, indicating a group velocity 
less than the vacuum speed of light c. In the region of anomalous dispersion near 
resonance the classical group delay is negative, indicating that the group velocity is 
greater than c. The steep slope of the index of refraction and the accompanying 
strong and rapidly varying absorption severely distort pulses that are not sufficiently 
monochromatic, and the standard physical interpretation of Auj/dk as the speed 
at which the peak of a modulation envelope travels is not valid for such pulses in 
this region. The "temporal-center-of-gravity" delay retains obvious physical meaning 
for all pulses, and in section [S] we develop quantum mechanical properties that are 
analogous to this delay for spontaneously emitted photons. This provides a framework 
for an unambiguous and causal interpretation of negative values of delays in both 
classical and quantum models. 

Group delays and "temporal-center-of-gravity" delays can be understood as a 
result of the transient oscillations in the medium [7]. These transient oscillations 
radiate fields that initially cancel the field of the incident wave, modifying the leading 
edge of transmitted pulses before they settle down to the steady-state fields given by 
equation Field transients must be included explicitly in the calculation of delays 
for pulses in the region of anomalous dispersion. We include a specific calculation of 
the "temporal-center-of-gravity" delay in section |S1 

3. Quantum Model 

The quantum mechanical system we consider is illustrated in figure [SJ and consists 
of three two-level atoms in a one-dimensional multimode optical cavity that extends 
from z = Q to z = L. The leftmost atom will initially be in the excited state, and 
will be the source of the field. The atom in the middle will scatter the radiation 
emitted by the source atom, and the atom on the right will serve as the detector. 
The cavity is assumed to be large in the sense that the length L is very much greater 
than the wavelength of the light emitted by the atoms. (The finite optical cavity does 
not contribute to the physical phenomena under investigation; it simply provides a 
convenient quantization volume for the field modes used in our calculations.) In this 
section we discuss the standard quantum optical Hamiltonian that we use, and we 
present the analytical solution for the time dependence of this system. 
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The zero- field resonance frequencies of the atoms are labeled where j — 1, 

2, or 3, and the positions of the atoms will be labeled Zj. In the results given below 
we will assume that the atoms are at positions zi = i/4, Z2 = L/2, and Z3 = 3L/4, 
as illustrated in figure [21 although the delay times we derive do not depend on the 
exact positions. The standing-wave field modes of the cavity are separated in angular 
frequency by the fundamental frequency 

Ac - 7r|. (13) 

This mode spacing may be small enough that many modes fall within the natural 
line-width of the atoms. 

For convenience we assume that the frequency of one of the modes corresponds 
exactly to the resonance frequency of atom 1, the emitting atom, and that the length 
of the cavity is such that it contains an even number of wavelengths of this mode. We 
label the frequency of this mode luq — mo Ac, where mp is an integer divisible by 4. 
(This assumption affects the details of some of our calculations, but not the existence 
of an analytic solution, nor any of our results concerning delay times. Classically, this 
assumption assures that optical path length between any pair of atoms is an integer 
number of half- wavelengths of the resonant radiation from atom 1.) The other mode 
frequencies will be enumerated from this mode so that 

ujm = (mo +m)Ac, (14) 

where m — 0, ±1, ±2, . . . Our large-cavity limit also assures that the frequencies of the 
three atoms are all very much greater than the fundamental frequency of the cavity, 
i.e., wj'^*^ 3> Ac. In this limit the atoms interact with a large number of modes, and 
not simply the resonant mode. 

As in the classical case we wish to study the effects of the detuning of the source 
field on the scattering of the radiation. We use the same symbol 6 as in the classical 
case to represent the detuning of the field, but in the quantum case the detuning is 
directly tied to the properties of the source and scattering atoms: 

S = 4"^'^ - uj^'^'K (15) 
The detector atom (atom 3) is assumed to have the same resonance frequency as the 

, . (at) (at) 

source atom, i.e., col = . 

We use as basis states the eigenstates of the atomic plus free-field Hamiltonian 

^^0 = -ff atoms + ^^ficld 
3 

= ^ huj^''^^aj + ^ hajmalr^am, (16) 

j=l m 

where fT| is the third component of the atomic pseudo-spin operator of atom j, and 
Om and al^ are the lowering and raising operators for the m^^ field mode. (We have 
re-zeroed the energy scale to remove zero-point energy of the field modes.) The basis 
states will be denoted as follows: 

• |e, g, g; 0) — Atom 1 excited, atoms 2 and 3 in ground state; no photons in field, 

• \g, e, g; 0) — Atom 2 excited, atoms 1 and 3 in ground state; no photons in field, 

• \g, g, e; 0) — Atom 3 excited, atoms 1 and 2 in ground state; no photons in field, 

• \9^ 5j 5; Im) — All atoms in ground state; one photon in field mode with frequency 
(mo -I- m)Ac. 
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We use the standard electric-dipole and rotating-wave approximations in the 
interaction Hamihonian |18[I19[I^ to give 

H -^atoms ^" -Afield ^ -^interaction 
3 

= Ho + J2J2^ {ajmamcrj + g*j„^ai,cJJ) , (17) 

j=l m 

where the strength of the couphng of the j"^ atom to the m^^ mode of the field is 
characterized by the constant gj^ , and cr^ and a~ act as raising and lowering operators 
for atom j . 

In our large-cavity limit, ijJj'^^^ S> A^, we can make the approximation that all 
modes that influence the dynamics of the system are near the atomic resonances, and 
the atom-field coupling constants are given by 

gjm — sin [(mo + jn)TTZj/L] . (18) 

In the previous equation Qj is a constant given by 

where dj is the dipole matrix element between the levels of atom j, and V is the 
effective volume of the cavity. 

Initially only the source atom will be excited and no photons will be present in 
the field, so that 

|V^(0)) = |e,5,5;0), (20) 
and we write the general state of the system as the linear combination 
IV'(i)) = ci{t)\e, g, g; 0) -I- C2(i)|g, e, g; 0) -f- C3{t)\g, g, e;0) +J2 bmit)\g, g, g; 1™). (21) 

m 

The Schrodinger equation yields the following set of coupled differential equations for 
the coefficients in equation (|21l) : 



, ,(at)„. 



9jmbyn (22) 



b,n = -i (Wm^m + 5l*mCi -|- 5^„C2 + g^^ca) . (23) 

We solve this set of equations with the Laplace transform technique used by Stey 
and Gibberd Laplace transforms have also been used to solve the Schrodinger 

equation in similar problems with two interacting atoms in three-dimensions |22l I23| . 
Because the Laplace transform technique is not new, and because we would like to 
focus on analogies with the classical model and physical interpretation, we leave the 
details of our solution to the appendix, and simply quote our results here. 

The general features of the solution giving the time-dependent atomic excitation 
amplitudes are illustrated in figure |3| The initially excited atom decays exponentially 
until t = 0.5L/c, the time at which scattered and reflected radiation flrst returns to 
the atom. The amplitudes to find the other atoms excited are identically zero until 
radiation first reaches them: the scattering atom first becomes excited at i = 0.25L/c 
and the detector atom is first excited at t = 0.5L/c. The three decay constants which 
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Figure 3. Magnitude of the amplitudes for the atoms to be in the excited state, 
starting from the state |i/'(0)) = \e,g,g;0). The decay rates of the atoms are 
71 = 4, 72 = 64, and 73 = 1024 (in the units of the figure) and the detuning is 
5 = 1.5672. 
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characterize the spontaneous emission rates of each of the atoms emerge naturaUy in 
terms of the parameters of the Hamiltonian as 



1] 



= i%r-- 



(24) 



The causal nature of the dynamics is evident in that all disturbances are propagated 
at the speed of light c via the quantum field. The abrupt changes in the complex 
amplitudes at intervals of 0.5L/c are a manifestation of the finite speed of light and 
the atomic spacing of 0.25L. 

The abrupt changes are manifested in our analytic solution for the complex 
amplitudes Cj{t) and bmit) as sums of terms with step functions that "turn on" at 
successively later intervals of 0.5L/c. In the following formulas giving these amplitudes 
we truncate the expressions so that only the first excitations of atoms 2 and 3 are 
included. We also note that the following equations are specific in some details to 
the atomic positions Zi used in this paper. The positions of the atoms enter via 
the coupling constants gjm, and changes in positions will result in changes in "turn- 
on" times and relative phases of terms. One manifestation of our large-cavity limit 
w ^ Ac is that the source atom decay at the free-space rates is clearly visible before 
interruption by reflected and scattered radiation, for example at times t < 0.5i/c in 
figure 13 The initial atomic decay rates given by equation (|24ll are not affected by the 
positions of the atoms; the effects of position dependence in the coupling constants 
compensate to keep the total initial decay rates independent of position. (Examples 
of the relationship between interrupted free-space decay and cavity-modified decay 
rates are given in j24j.) Later in this paper we will consider decay times that are 
much shorter than L/c so that the effects of all initial excitations are complete before 
reflected or mutiply scattered radition excites the atoms. The conclusions of this paper 
concerning delay times are unaffected by any details of position dependence. 

The complete time dependence of the system is given by the following set of 
amplitudes (in which the time is scaled in units of L/c): 



ci{t) = exp 



2 



e 



exp 



1 

t 

2 



71 



(71 

1 

2 



72 + 12(5)2 

■(7i- 272 +i2<5)(7i-72 + i2(5) 



-I-72 exp 



(11 
V 2 



C2{t) = e(t 



71 — 72 + i' 



2S 



iS] t-^ 



71 



(25) 



exp 



exp 



csit) = Q[t-^] V7173 <! exp 



^ll(t-l 
2 V 2 



brn{t) = 



— exp 
+ exp 

i2gim 
(71 — i2m7r) 



2 V 2 



exp 



(71 + 125) 



(71 - 73) (71 - 72 + 12(5) 

72 

(71 - 72 + 125) (-72 + 73 + 125) 
(73 + 125) 



(-71 + 73)(-72 + 73 + 125) 



(-yi) - exp {-imnt) + Q U - ^ ) i2g2mV7i72 



(27) 
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exp [-l(t-i)] 



exp 



(71 — i2m7r)(7i — 72 + i2(5) (71 — i2m7r)[72 — i2(5 + m7r) 



exp 



[72 — i2((5 + m7r)] (71 — 72 + i2S) 



(28) 



(We note that if we had made a different assumption about the size of resonant 
wavelength relative to the length of the cavity in our discussion immediately preceding 
equation (|14|l . then the complex phases associated with terms in the equations above 
would be different.) 

In the following sections we will focus on two quantities: C3{t), the amplitude to 
find the detector atom excited, and the expectation value of the square of the 
electric field operator, which is proportional to the field intensity. (The expectation 
value of the field operator itself is zero for any state with the form of equation l|21|) .) 
In our investigation of c^it) we will consider only the displayed term in equation (|27|l 
describing the initial excitation of the detector atom. Similarly, we will investigate 

in regions to the right of the scattering atom, and at times that exclude multiple 
scattering effects. 

It is useful to rewrite C3{t) as the sum of two pieces: the amplitude c'^{t) for atom 
3 to be excited in the absence of the scattering atom (or, equivalently, when 72 =0), 
and c|(t), the amplitude that is attributable to scattering. The total amplitude is 
thus 



csit)^cl{t)+clit). 
Setting 72 = in equation l(77|l gives 

2 



lit)^eit 



exp 



t 



71 - 73 

and subtracting this from equation H27|) gives 
1\ f exp[-f(i- 



exp 



73 

' 2 



72\/7i73 



exp 



0] 



exp[-(f -i<5) {t 



(29) 



(30) 



(71 " 73) (71 - 72 + i2(5) (71 - 72 + i2(5)(-72 + 73 + i2i5) 



{t - ^)] 



(-71 + 73) (-72 + 73 + i2(5) 



(31) 



As an alternative to finding the time-dependence of the excitation amplitude for 
the detector atom, we can characterize the transmitted field itself without recourse to 
the details of the detector. Standard photodetection theory [21] suggests calculation of 
the expectation value (£^~-'(z,t)£'+^(z,t)), where £'+^(z,t) and £^~\z,t) correspond 
to the decomposition of the interaction representation field operator into positive and 
negative frequency parts. This is equivalent to the calculation of the expectation value 
of the normally ordered intensity operator (: £^ :). 

Using the electric field operator in the form given in reference ^HIj we write the 
expectation value of the square of the field as 

' \ 1 /2 ^ 



^oV J 



{an 



TT Z 

(mo + m) — 



\m)- 



(32) 



In the limit considered in this paper we can replace the frequencies lo„i under the 
radical with the constant uj^^^^ . After expanding the state vector as in equation H21f) , 
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normally ordering the operators, and evaluating the sums, the expectation value can 
be written in terms of the amplitudes b„i{t) to find the photon in the various cavity 
modes: 

2 

^5m(t)sin (mo + m)— . (33) 




Evaluation of this expression gives a space- and time-dependent representation of the 
localization of the energy of the photon |13l III] . 

The expression for (: £^ :) in equation 1)33(1 is the square of a complex number 
that is analogous to the classical complex analytic signal. We label this quantity Eq.m, 
i.e., 




(t) sin (mo + m)- 



L 



(34) 



The overall phase of this quantity is clearly arbitrary; in what follows we retain the 
phase that comes from a direct evaluation this equation. 

We note that quantity £q.ni. is very closely related to what has been identified 
as "the 'electric field' associated with [a] single photon state" by Scully and Zubairy 
|25|. They consider states which are the product of separable atomic and field states. 
They conclude that for single photon states 1-07) the quantity {0\&^^\'ip^) "can be 
interpreted as a kind of a wave function for a photon." They also demonstrate the 
classical nature of this quantity for field states produced by spontaneous emission. 
The present paper extends this kind of analysis to define an analog to the classical 
analytic signal even in cases in which field and atomic variables are entangled, and 
it extends the quantum-classical correspondence to fields that include scattering from 
two-level atoms. 

For ease of comparison with previous results for the detector atom, we fix 
z = 3L/4 in what follows. With no scattering atom present we find (see Appendix) 

(at) \ 

null Ti"" 



2eoVAc 



exp 



(35) 



The energy density passing the point z — ZL/A exhibits an abrupt turn-on (because 
of the initial conditions we have chosen) followed by exponential decay ^1 114L . 
With a scattering atom present at z = L/2 we find 

1/2 



-ie t 




exp 



71 
' 2 



72 



exp 



[26 - i(7i 
1^ 



72 j 



(36) 



Using only the first exponential term leads to the previous result with no scattering 
atom present; the effect of the scattering is contained in the remaining terms. 
Calculation of the expectation value of the square of the field will exhibit interference 
between the terms of this expression, leading to effects similar to those noted in 
reference [T^. 
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4. Comparison of classical and quantum mechanical scattering 



In order to compare the quantum mechanical "pulses" derived in section |3 to classical 
analogs it is necessary to extend the treatment reviewed in section El to include 
the effects of transients. This can be done in the frequency domain using Fourier 
techniques, or in the time domain by explicitly including transients in the solution of 
the equation of motion of the charged oscillators, and calculating the fields re-radiated 
by the transients. The quantum mechanical results of section |31 have been derived in 
the time domain, and to emphasize the analogy between the classical and quantum 
cases we review the classical time domain calculation. 

The quantum pulses of equation H35|l that are incident on the scattering atom 
have the classical analog 



^ci. = U - 77 C-exp 



— + 10^1 



1 

t 

2 



(37) 



where C is a constant. Using the field of equation (|37|l as the driving field in the 
equation of motion of a driven charged oscillator with damping constant 72 and 
resonance frequency uj2, it is straightforward to find the motion of the oscillator. 
For the initial conditions a::(0) = and i(0) — 0, and assuming that uj2 2> 72, we have 



x{t) ~ e 



1 

'--2 



qC/m 



{ujj -ujI+ iwi(7i - 72) + 7i/4 - 7172/2) 



exp 



/7i 

vy 



ICiJl 



1 

'--2 



exp 



1 

*-2 



(38) 



Calculating the re-radiated field from the motion of the charges as in equation I^J, 
and making the approximations (5 = o^i — 0^2 ^ l^j and 7^ ^ w, gives the scattered 
field as 

- ) exp (-iwit) 



-ie t 



[2<5-i(7i-72)] 



exp 



71 
' 2 



exp 



t 



(39) 



We note that the functional form of this equation giving the classical scattered field 
is identical to the portion of equation l|3t)|) describing the quantum scattered field. In 
addition, we note that this functional form is also refiected in c^lt). In the limit of 
rapid detector response, i.e., 73 3> 71,72, this amplitude is 



C3(0 



-e t 



exp 



exp 



71 
' 2 



t 



71 
' 2 

1 
2 



1 

^ 2 
exp 



72 



[2(5 - i(7i - 72 
'72 .A 1 



(40) 



and the contribution due to scattering is again the same. (Our use of a "detector" 
which is an atom identical to the "source," except with a more rapid decay rate, is not 
meant to correspond to any real experiment. It is an idealization meant to illustrate 
the effects of the quantum field on a simple system with a fast response.) 

Thus we see the exact functional form of the classical field reflected in two 
quantum mechanical quantities: the probability amplitude of the detector atom C3(t), 
and in our quantum field "amplitude" £q.m.- The classical and quantum expressions all 
include a "steady-state" term (which drops off at the slow decay rate of source atoms. 
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7i/2). In the limit 71 <C 72 the radiation emitted by the source has a very narrow 
line-width compared to the scattering atom, and this steady-state term simplifies to 
the form of equation {TJ given by the simple classical theory of section |2 

The expression for the quantum field "amplitude" Sq.m. contains a time-dependent 
phase that is identical to the phase of the classical field, and the probability amplitude 
Csit) contains terms due to incident and scattered fields with relative phases that 
match classical expectations. The concept of the phase of a quantum field has been 
the subject of investigation from the early days of quantum mechanics jSSlEZl- We 
note, however, that most previous work on the phase of a quantum field has focused 
on defining a meaningful and mathematically well-behaved phase operator for linear 
combinations of multiply occupied states of a single-mode field; in this paper we see 
manifestations of classical phases in a multimode field which contains a single photon. 

In the case of classical scattering considered in section |21 the effect of a single 
scattering event is considered to be small, and exponential phase rotation and 
attenuation of the transmitted field are the result of many scattering events. In our 
quantum model the magnitude of the scattering is determined by Q2 (or equivalently 
72) which characterizes the coupling of atom 2 to the field. In our one-dimensional 
model the coupling to the incident field and the decay rate of atom 2 are both 
completely determined by the single parameter 0,2 , which means that it is not possible 
to make the effect of the scattering small without simultaneously making the line- 
width of atom 2 very narrow. In a fully three-dimensional model the decay rate of 
atom 2 would be the result of the atom's coupling to many more modes, and not just 
those containing the incident field. This scattering into other modes would reduce the 
scattering in the forward direction (the direction of the detector) from the amount 
predicted in our simple model. Of course in a fully-three dimensional model there 
are other effects also: the fraction of source radiation emitted in the direction of 
the scatterer (and detector) would be reduced, and the scatterer would be driven by a 
weaker field and would emit less total radiation. The net effect is that the excitation of 
the detector due to the scattered field will be reduced relative to the direct excitation. 
If the relative effect of forward scattering is reduced by a factor /, then a more realistic 
expression for the excitation of the detector atom is 

cs{t)=U{c°{t) + fcl{t)), (41) 

where Af gives an overall reduction in the excitation. We have used the same symbol 
/ here that we used earlier for the dimensionless parameter which characterizes the 
magnitude of classical scattering in section[2| The quantum superposition of equation 
(|41|l is analogous to the classical field superposition of equations © and . 

5. Temporal-center-of-gravity delay 

The group delay of a classical pulse has a clear interpretation for quasi-monochromatic 
pulses: it is the delay in the arrival of the peak of a pulse compared to the time expected 
for propagation through a vacuum. The pulses investigated in this paper have sharp 
leading edges, and this lack of a smooth modulation envelope means that the results 
of simple classical theory for quasi-monochromatic pulses should not be expected to 
be a sufficient guide to full understanding. In this section we investigate "temporal- 
center-of-gravity" delays in several classical and quantum mechanical quantities. 

The first delay we investigate is derived from C3{t), the amplitude for the detector 
atom to be excited. As we have argued previously, this amplitude will reflect the 
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Figure 4. Probability for the detector atom to be excited as a function of 
time. Tlie top figure is for a detuning 5 in the region of normal dispersion, 
and the bottom is in the region of anomalous dispersion. The effect of the 
"pulse re-shaping" is different in the two graphs. In the top graph probability 
is preferentially removed from earlier times, while the converse is true in the 
lower graph. The decay rates of the atoms are 71 = 1, 72 = 16, and 73 = 256 (in 
the units of the figure). 

strength of the incident field in the hmit that the response time of this atom is very 
small compared with other time scales, i.e., 73 71,72- The effect of the scattering 
on this amplitude is evident in figure 0] in which |c3(i)p is plotted for two values of 
the detuning 5, and also for the case in which no scattering atom is present. (For 
clarity in this figure we have only included the terms describing the initial "turn-on" 
of excitation; we have not included effects due to reflection and multiple scattering.) 
For all detunings, C3{t) = for all times earlier than t = 0.5L/c, as is expected; all 
effects on the detector atom occur at times that preserve causality. The qualitative 
shapes of the detector response depend critically on the detuning 6. For detunings 
in the region of normal dispersion the response is relatively reduced at early times 
compared to later times, resulting in a qualitative delay in detection of the photon. 
For detunings in the region of anomalous dispersion (—72/2 < S < 72/2) the converse 
is true. 




5 = 72/4 
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We quantify these ideas by identifying an effective arrival time of the photon with 
the temporal center of gravity of the probability that the detector atom is excited, 
i.e., 

_ jt\cm'dt 

tarrival- J \^^^t)\2 ' ^^^i 

(In evaluating the integrals in this equation, we use only the term in the series of 
equation (|27|l that "turns on" a,t t — 1/2, and assume that the decay rates and 
distances are such that the effect of multiple scattering is negligible. These decay rates 
are very much faster than those displayed in the figures; the rates in the figures were 
chosen so that the the atomic excitation dynamics and the causal delays occurred on 
the same time scale, and could be illustrated on the same graphs.) The delay imposed 
by the medium is then just the difference in the arrival times calculated with and 
without a scattering atom present, 

_ Jt\cm'dt /t|cg(Opdt 

J\c,it)\'dt JlcUtWdf ^^'^ 

To explore the effect of weak scattering we rewrite C3{t) in the form of equation 141|) . 
and assume that / ^ 1. Our quantum mechanical delay becomes, to first order in /, 
^ /^|c°ft) + /c|ft)pd^ Jt\c°,{t)\^dt 



2/ 



JtRe [cl{t)cl{tY] dt /t|c§W|2di/Re [c§(t)c|(t) *] dU ^ 



Atc3-4/7. ;;;, , % . (45) 



!\cl{t)?dt (/|c§(i)pdt)' 

It is straightforward to evaluate the integrals in equation H44|l using the expressions 
for c%{t) and c\{t) from equations lp!n|l and lp?T|l . After taking the limit 73 ^ oo and 
then letting 71 ^ we arrive at the following expression for the quantum time delay: 

(4<5^-7l) 

Comparing this to equation (|12l) shows that "temporal-center-of-gravity" delay time 
for this specific pulse is identical in functional form to the classical group delay. The 
magnitude of the "temporal-center-of-gravity" delay is, however, twice that given by 
the group delay. (The "temporal-center-of-gravity" delay for classical pulses with the 
form given in equation (|37|l is also twice as large as the group delay.) 

The delay in the arrival time that we have defined is the result of the reshaping 
of the "pulse" of excitation of the detector atom. The effect of scattering in the 
region of anomalous dispersion is to reduce preferentially the probability that the 
detector atom will be excited at large times; a short spike of probability at early times 
remains, shifting the "center-of-gravity" of the excitation to earlier times. Despite the 
appearance of negative delays, no excitation of the detector atom occurs at times 
t < 0.5L/c for any values of the parameters in our model. Although we have 
not demonstrated it explicitly in this work, we are confident that delays of quasi- 
monochromatic quantum pulses can be explained in the same manner. 

Because £ci., i^^q.m., and 03(1) all have the same functional form (in the large 73 
limit) it is easy to see that equivalent delays can be derived from the classical field 
using equation Q, or from the quantum field using the analogous equation 

_ f jt{:£iz = 3L/ir :)dt \ / /t(:g(z = 3L/4)^ :) dA 
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6. Conclusion 

We have considered simple classical and quantum models of propagation of light 
through dispersive media, and we have identified quantum mechanical quantities that 
have exact analogs in the classical scattered field. Scattering induces delays in the 
quantum model which are identical to the delays of a classical field. Because we use 
standard field-theoretic techniques of quantum optics to account for the creation and 
absorption of photons in our model, we have a well-defined initial condition in which 
all of the energy is localized at the position of the initially excited atom. This enables 
us to demonstrate that no effects propagate faster than the vacuum speed of light c, in 
spite of the appearance of negative delays that seemingly correspond to superluminal 
group velocities. In a companion paper we apply similar techniques to study delays 
of photons in media exhibiting electromagnetically induced transparency |28| . 

Appendix A. Solution using Laplace Transforms 

Taking the Laplace transform of the coupled differential equations 122|I - H23() gives the 
coupled algebraic equations 

i(sci(s) - 1) = ^5im6„(s), (A.l) 

711 

isc2(s) = -6c2is) + ^ g2mbm{s), (A. 2) 

m 

isc3(s) = ^53„6„i(s), (A. 3) 

m 

isbjn{s) = mAc6m(s) +^Cj(s)g*„. (A. 4) 

i 

The idea behind our solution is straightforward: solve these algebraic equations 
for the quantities Cj{s) and bm{s) and then perform an inverse Laplace transform 
to recover the time dependence of Cj{t) and hm{t). The details of carrying out 
such calculations are quite involved, and were completed with the aid of a computer 
algebra system. + In this appendix we outline our approach and present some of our 
intermediate results. 

We begin by solving equation (jA.4|l for and substitute the result in the 

first three equations, giving 

sci - 1 = -iAc (/iici -I- /12C2 -I- /13C3) , (A. 5) 

(s - i6)c2 = -iAc (/21C1 + /22C2 + /23C3) , (A. 6) 

SC3 = -iAc (/31C1 + /32C2 + /33C3) , (A. 7) 

in which we have defined the dimensionless sums 

r 1 \ ^ 9lm9nm I \ Q \ 

*" m Ac 

In the limit in which the atomic resonance frequencies are very much greater 

(at) 

than the fundamental frequency of the cavity, i.e., utj 3> Ac, these sums may be 

+ Mathematica notebooks used to perform the calculations are available from the authors. 
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approximated by extending the range for m from — oo to +00, in which case the 
sums have relatively simple representations in terms of trigonometric functions. An 
explicit example of one of these sums (for atoms at positions zi = L/A, Z2 = ^^/2, and 
Z3 = 3L/4) is 



8a: 



cot 



-\ -. 1 + 2 cot 



- 4^-2 



cot 



IS 

a: 



71 
4A, 



sinhjl^ 
cosh cosh 



After solving equations ljA.5p - (|A.7p for the quantities Cj{s) in terms of the sums 
/jm, we rewrite the hyperbolic trigonometric functions resulting from the sums in 
terms of exponential functions; we also let c/ L — 1 at this point in the calculation. 
We then expand the resulting expressions in powers of exp(— s/4), and the time 
dependence of the system is recovered by a term-by-term inverse Laplace transform 
of the expansion. The step function turn-on of the resulting time dependence arises 
because of the factors exp(— ns/4) in the expansion. The lowest order terms in our 
expansions of the Laplace transforms are given here: 

2 4exp(-s/2)7i(s 72 - i(5) 



ci(s) = 

C3(s) = 
bm{s) = 



2s + 7i (2s + 7i)^(2s + 72 -i2(5) 
2cxp(-s/4)^7i72 



(2s + 7i)(2s + 72 -125) 
4 exp(-s/2) ^7173(5 



i<5) 



(2s + 71) (2s + 73) (2s 72 - i2<5) 



2i 



cxp(-s/4)^7i72g2„ 
(2s 72 - 12(5) 



9ln 



q.m. of equation (|36|l we repeat the process with the detector atom 



(s + im7r)(2s -|- 71 

The inverse Laplace transform of these expressions gives equations H25|l -p8 |) 

To calculate 
at z = 3i/4 removed, yielding 

ci(s) = 

C2(s) = 

^m(s) = J. TTTT, \ T 1 9lm - 92r 



(A.IO) 
(A.ll) 

(A.12) 
(A.13) 



2s- 



71 



2exp(-s/4)^7i72 



(2s- 



-71) (2s 
2 



72 



(is-TOAc)(2s-f 71) 



m) 

9ln 



^7i72exp(-s/4) 
■ (2s- 



(A.14) 
(A.15) 

-(A: 16) 



Using the same techniques as above, the Laplace transform of Sq.m 

/ m37r\ 



72 - i2<5) 

can be written as 



eoV 



Xs) sin 



V 4 



^2?lu;^'''^7iexp(-s/2) 
TrAceoV {2 s + ji) 



72 



(2s -t- 72 - i2<5) 



(A.17) 



The inverse Laplace transform of equation (jA.17|l gives equation 
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